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Метод адаптивной рекуррентной фильтрации 
измерений параметров космических объектов 


В работе рассмотрен синтез адаптивного рекуррентного фильтра, который предназначен для 
статистической обработки измерений параметров движения космических объектов на этапе 
орбитального маневра в измерительных информационных системах. В основе разработки лежит 
принцип разделения. В фильтре учтена особенность модели движения космических объектов, что 
позволило эффективно осуществлять идентификацию параметров траектории на этапе орбитального 
маневра, а также существенно снизить динамическую ошибку сопровождения на данном этапе. 


Постановка проблемы 


Известно, что для определения и уточнения параметров орбит космических 
объектов чаще всего используются измерительные информационные системы, которые 
применяют для измерения текущих координат радиолокационные станции или 
лазерные локаторы. При этом для обработки поступающей информации целесообразно 
использовать рекуррентные фильтры, поскольку количество траекторных измерений 
обычно бывает достаточно большим (сотни-тысячи измерений) и хранить их в памяти 
ЭВМ нерационально. Наиболее известным и широко используемым является фильтр 
Калмана [1]. Известны ряд модификаций этого фильтра, которые позволяют устойчиво 
сопровождать космические объекты на детерминированных участках траекторий. 
Однако на участке орбитального маневра космического объекта подобные фильтры 
приводят к появлению недопустимо больших динамических ошибок сопровождения, 
которые приводят к сбросу с сопровождения таких объектов. 


Анализ публикаций 


В работе [2] предложен метод регуляризации рекуррентного фильтра, состоящий в 
добавлении к матрице ошибок измерений дополнительного члена, что позволяет на 
участке слабого маневра осуществлять сопровождение космических объектов, однако 
при этом невозможно осуществить оценку истинных параметров интенсивности 
маневра, а также присутствует достаточно болышая динамическая ошибка сопровожде- 
ния, что может привести к сбросу с сопровождения объекта. 

Для решения указанной проблемы применяются фильтры с обнаружением 
маневра, оценкой его интенсивности и переключением на рекуррентный фильтр с 
соответствующей моделью движения. Однако такой подход часто требует больших 
затрат вычислительных ресурсов и не позволяет существенно снизить динамическую 
ошибку сопровождения на переходных участках. 

Другой подход состоит в применении адаптивных рекуррентных фильтров. 
Предложенный в работе [3] адаптивный фильтр позволяет эффективно решать задачу 
сопровождения объектов, которые осуществляют маневр относительно линейной траек- 
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тории, при этом динамическая ошибка сопровождения таких объектов снижается 
в 2—3 раза по сравнению с фильтрами с обнаружением маневра. Однако примене- 
ние такого адаптивного фильтра для решения задач обработки данных траекторных 
измерений космических объектов не представляется возможным, поскольку в рассма- 
триваемом фильтре базовой является линейная модель движения, а модель движения 
космических объектов на достаточно длительных интервалах сопровождения не является 
таковой. В результате фильтр постоянно функционирует в режиме сопровождения 
маневрирующего объекта, что приводит к появлению дополнительной динамической 
ошибки. 


Цель статьи 


Осуществить синтез адаптивного фильтра, в котором учитывается особенность 
модели движения космических объектов, для обеспечения возможности эффективной 
оценки параметров орбитального маневра и устойчивого сопровождения измери- 
тельными информационными системами на таких участках орбиты. 


Формализация задачи 


Для успешного сопровождения измерительными информационными системами 
космических объектов на участке орбитального маневра необходимо использовать 
модель движения, в которой присутствуют дополнительные члены, описывающие 
такой маневр. 

В статье используются обозначения переменных, которые введены в [3] и других 
работах С.З. Кузьмина, поскольку им наиболее полно рассмотрены практические вопросы 
рекуррентной фильтрации в измерительных системах. 

Модель движения космического объекта на участке орбитального маневра 
может быть представлена в виде 


о, == Ф о г Гб +С 77, 2 (1) 


пп 


где Фо, - уравнение базовой невозмущенной траектории (для космических объектов 
целесообразно использовать полином второй степени), ®„ — [-мерный вектор возму- 
щений на участке орбитального маневра космического объекта, 77, — р-мерный вектор 
возмущений, вызванных влиянием внешней среды, Г, и С, являются известными мат- 
рицами, позволяющими перевести элементы векторов #„ и 1, в единую систему коор- 


динат, используемую в модели движения. 

Следует отметить, что измерение параметров траекторий космических объектов 
измерительными информационными системами (например, радиолокационными 
системами) обычно осуществляется в сферических системах координат [4], однако в 
них даже невозмущенная траектория движения описывается уравнениями, содержащими 
производные порядка выше двух, что чрезвычайно усложняет рекуррентный фильтр, 
увеличивает его размерность и делает его нереализуемым в реальном масштабе времени. 

Для упрощения фильтра в работе [3] предложено фильтрацию осуществлять в 
топоцентрической системе координат [4], а для поддержания сопровождения цели 
осуществлять пересчет экстраполированных координат в сферическую систему 
радиолокационной станции. 
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Таким образом, синтезируемый фильтр должен обеспечить рекуррентную фильтра- 
цию поступающих измерений параметров космического объекта, уравнения движения 
которого описываются соотношением (1), а также обеспечить эффективное сопровож- 
дение объекта на участке орбитального маневра, когда значения элементов вектора 
2 „ априорно неизвестны. 


Синтез адаптивного рекуррентного фильтра 


На этапе орбитального маневра возмущение траектории (1) космического 
объекта может быть представлено случайным процессом [3], среднее значение 
которого может принимать значения из фиксированного множества состояний в 


диапазоне [-2к›...›0,....+8 „пик |. Переходы из состояния { в состояние / совер- 
шаются с вероятностью л,, время пребывания в состоянии 1 является случайной 
величиной с плотностью вероятности @(г). Такой процесс является полумарковским 


случайным процессом [3]. 
В работе [5] сделан вывод об эффективности принципа разделения для построения 
адаптивных систем. Целесообразно его применить для построения адаптивного фильтра. 
При оценивании вектора параметров движения космического объекта 0, 


предположим, что известна апостериорная плотность вероятности @ (о, | {%},) вектора 


х, по данным п -мерной последовательности измерений {У}. Тогда при квадратичной 


функции потерь апостериорное математическое ожидание (оценка) может быть вычис- 
лена по формуле Байеса 


Жо | 0,6 (ь, {У}, а, . (2) 


Здесь (9) — пространство возможных значений оцениваемого параметра. 


В случае присутствия возмущающего воздействия апостериорная плотность 
вероятности может быть представлена в виде [6] 


о (>.1(У,)= ов. У, У, 48... (3) 


Здесь (=) — область возможных значений интенсивностей орбитального 


маневра. 
Осуществляя подстановку (3) в (2), получим соотношение 


6, = [ь, [ ©(6, |5 „(У (2 „„ {У}, ) 4 „Чо, (4) 
@ @) 


м 


После преобразования выражения (4) имеем 


6, = |6, (8„)0(5 „| {У}. (5) 
(=) 
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В работах [3], [6] предложен оригинальный подход к вычислению значений 


плотности вероятности «(5 | {У}, ) ‚ которые необходимы для получения значения 


0, ‚ то есть решения задачи адаптивной фильтрации. 
Поскольку параметр ®„ принимает только фиксированные значения, а также 


учитывая свойства определенного интеграла [7], выражение (5) преобразуется к 
виду 


[2 
0, № о 0, [Е (5 | {У}. ) 2 (6) 
р 
где #„„ — одно из ] (1= —К,— -(&-1)....,0,1,...,К-ЬА) фиксированных значений 
параметра $; 
Р(. {У}, — полученная по результатам и измерений {У} апостериорная 


вероятность события 2, =. 


Применим правило Байеса [8] для вычисления Р(. {У}. ) ‚ тогда 


Р(з {У}, .) (У, 18.) 


ТТ о (7) 
. | Хо Р(в {У} (У, 18 ь-о) 


1% 


Здесь Р(. {У} — рассчитанная по результатам п-1 измерений априорная 
вероятность параметра маневра #„ на п-ом шаге фильтрации. Для вычисления 


этой вероятности, учитывая полумарковость процесса маневра, применимо соот- 
ношение [3] 


Р(е„.|{У}, = У л,Р(в и: (8) 


1=—2К 


где л, =Р (=. — 8 |8 = 5.) — условная вероятность перехода процесса маневра 
из состояния # на (и-1)-м шаге в состояние / на и-ом шаге. 
Условная плотность вероятности о(У, | [а в измеренного значения век- 


тора параметров У, ‚ если параметр орбитального маневра на (п -И-м шаге имел 
значение $ „,, обычно аппроксимируется нормальным распределением [3], [6] с мате- 


матическим ожиданием У„, и дисперсией у’ и вычисляется по формулам 


У =Н, [Ф0,+Г, 5. |, (9) 


эп 
ЕН, „НТ+0?. (10) 
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Здесь Н,, Ф,/ - известные матрицы наблюдения и экстраполяции фильтра 
Калмана. Тогда учитывая соотношения (7) — (10) и обозначив Р(. {У}, =» 
имеем выражение для апостериорной плотности вероятности 

^ ^ 2 
Е 
п; (= п- )схр ОЕ ЗЕ Е: "=. 
и И) м(п-—1) п-1 2а 
‚РА (11) 


В (У, -У, 
У ПР У, )ер Р. 


УЕ 


2 
п 


Из анализа полученного соотношения (6) может быть сделан вывод о том, что 
практически задача оценивания искомого вектора 0, сводится к взвешенному 


усреднению оценок 0, [Е ) ‚ которые могут быть получены любым известным спо- 


собом при фиксированных значениях © 


мт * 
Для получения оценок 0, (*„) целесообразно применить фильтры Калмана, 


поскольку они традиционно используются в измерительных информационных 
системах и хорошо себя зарекомендовали. Тогда синтезированный адаптивный 


фильтр представляет собой набор фильтров Калмана для получения оценок о, (и ) у 
которые на выходе суммируются с весовыми коэффициентами Р, для нахождения 


искомой оценки о,. Структурная схема полученного адаптивного фильтра показана 


на рис. 1. 
Набор операций, которые выполняются в фильтре Калмана, известен [1], [3], [6] и 
может быть представлен в виде следующих выражений 


о 58 Ф, о 1} 
ее = ФР, .Ф,; 
К, =Ч4„Н! (Н.Н! +В, ); (12) 


5, =5„+К, (У,-Н,5 


Е = та -К НЧ. 


т. 


Из анализа соотношений (12) можно сделать вывод, что значения коэффициен- 
тов усиления фильтров К, , а также ковариационных матриц ошибок оценок пара- 


метров Ч, не зависят от конкретного значения элементов вектора возможных 
значений интенсивностей орбитального маневра $,„„. Таким образом, эти параметры 


могут рассчитываться в отдельном блоке и использоваться сразу во всех из 2+1 
фильтрах Калмана. Учитывая данное обстоятельство, можно сделать вывод, что 
вычислительные затраты, которые необходимы для реализации рассматриваемого 
адаптивного фильтра существенно снижаются, что является важным фактором для 
систем обработки информации реального времени. 

Тогда структурная схема адаптивного фильтра [3] может быть представлена в 
следующем виде (рис. 2). 
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ГлАы 
Е ЗВАЖЫЫ ` 
> , су 
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Гы 
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Рисунок 1 — Общая структура рекуррентного адаптивного фильтра 


пу = +, 
] в а 
=> к + > Блок Р ил | >| Блок А ы 
в ТА памяти ‚| памяти Г 5 ат 
п п 
у 
> 
_| (ии) ГХН 8 
ИО 
У, у 
Ао . Т.Н 
_> © (+) ы 
и ^ 
|” = о, 
5, (Е) Хх а › 
Гл 
у 
РУ Е НИЧ ) 
© о Ре . 
г > р 
ы ОИК 2 
О а Блок 
о ен $, (8) хе памяти 
А я 
5, ты о 


Рисунок 2 — Структурная схема адаптивного фильтра 
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Рассмотрим содержание основных операций, которые реализуются в адаптивном 
фильтре, с учетом ранее указанных особенностей обработки информации о параметрах 
космического объекта на участке орбитального маневра. 


Пример реализации адаптивного фильтра 


Следует отметить, что измерения координат космических объектов радиолока- 
ционными станциями обычно производятся во второй сферической системе координат [4] 
К, =, у. Однако вследствие особенностей такой системы ускорение может достигать 
существенно больших величин, чем в прямоугольных системах координат. 

К примеру, при скорости полета И = 8 км/с на дальности А = 500 км при №=0 
значение радиального ускорения составит Ё = 128 м/с, в то время как в прямоугольных 
инерциальных системах координат максимальное значение ускорения не превышает 
величины ускорения свободного падения ©, = 10 м/с”. Это приведет к тому, что фильтр 
будет работать постоянно в режиме маневра, кроме того, потребуется больше блоков, 
которые реализуют фильтры Калмана, которые настроены на определенную интенсив- 
ность маневра. Следовательно, обработку информации в адаптивном фильтре целесо- 
образно производить в прямоугольной системе координат. 

Применим для этого топоцентрическую систему координат [4], как непосред- 
ственно связанную с местом размещения измерительной информационной системы. 
Начало ее находится на поверхности Земли (Землей понимается ее шаровая модель). Ось 
‘у направлена в зенит по нормали к поверхности, а оси х и 2 лежат в горизонтальной 


плоскости. Направление оси 2 задается азимутом А, (А, — угол, отсчитываемый от 


направления на север по часовой стрелке до оси 2). Ось х дополняет систему до правой. 
Топоцентрические прямоугольные координаты выражаются через сферические 
К, =, у при помощи соотношений 


х = Аз &с03$у; у= АзшЕзЗШУ; 2=Ас035. (13) 


Здесь & (0<=<л) - угол между осью 2 и линией, соединяющей начало коорди- 


нат с точкой а(К, =, У). Угол у (-л<у<л) представляет собой линейный угол 


двугранного угла, образованного горизонтальной плоскостью и плоскостью, прохо- 
дящей через ось 2 и точку а. 


Будем считать, что на (и- 1) -м шаге получены вектор оценок параметров траекто- 


2 т 2 р Е В т 
рии космического объекта 0, о О М ) ‚ апостери- 


27 п-—1 
орные вероятности возможных значений интенсивностей маневра Р(®., {У} ) 
( 1=-=&,- (к — 1) ‚..., 0,15... -ЦЬ К) ‚ ковариационная матрица случайной составляющей 


и шумов управления Ч „,„_), а также ковариационная матрица ошибок оценок 


9 Уре ^” Иво Ив) 
к) >20 ^° 30 Иа 
= : : 5% : - . (14) 
Уч У ^ Из Из) 
-) и-) ^° и и) 
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При этом следует учесть, что элементы Уз, (1=1....,9) являются дисперсиями 


ошибок оценок соответствующих параметров вектора [9 а остальные элементы 


п-1? 
представляют собой соответствующие ковариации. 
Экстраполяция параметров траектории ДЛЯ каждого ИЗ ВОЗМОЖНЫХ 


о (1 =-,- (& 1,...,0,1,...&-Ь®) осуществляется по формулам 


Вы (=„)= Фо, 1 +Ф и. 2 (15) 
где матрицы Ф,, Ф,„ — матрицы экстраполяции векторов состояния и маневра 
соответственно. При этом 

Ф 0 0 
Ф= о Ф 0, (16) 
ооф 
0 че / 2 
где Ф.=| 0 1 г, | матрица экстраполяции части вектора состояния, соответст- 
оо 1 


вующей определенной координате и ее двум производным. 
Матрица Ф„ имеет вид 


Ф„ 0 0 0 т. / 2 
Ф„ = 0 Ф.:‹ 0 ? Ф„ = 0 0 т, 2 ( 1 и) 
0 оф оо 


гдет тЫ интервал экстраполяции, а расширенный вектор возможных интенсивностей 
И 


Расчет ковариационной матрицы ошибок экстраполяции осуществляется по 
формуле 


маневра 5 = (0, 0х, 
Ф,=ФЧ, Ф". о 


Вычисление вектора коэффициентов усиления фильтра в соответствующем 
блоке производится в соответствии с выражением 


К, =Ф,Н! (Н.Ф „Н! +В, , (19) 
0о00000о00 
где матрица наблюдения Н, =| 0 0 0 1 0 0 00 0, ковариационная матрица 
000000100 
о, 0 
ошибок измерения В =| 0 с, имеет диагональный вид, поскольку в 
0 | 
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соответствии с рекомендациями, которые изложены в работе [3], целесообразно 
перейти к упрощенному методу фильтрации и не учитывать ковариации, которые 
возникают при пересчете матрицы ошибок измерения из сферической системы 
координат в прямоугольную. При этом потеря точности не превышает 5 — 15 процентов, 
но существенно снижаются вычислительные затраты. 

Ковариационная матрица ошибок оценок параметров определяется по формуле 


ен. = ее ИО кН,» А (20) 


Оценки фильтруемых параметров для каждого из возможных значений 
интенсивностей маневра 2, вычисляются по формуле 


5, =0,„(8„)+К, (У,-Н,5, (8), /=-№...0,..5К. (21) 


Апостериорные вероятности (веса возможных значений интенсивностей 
маневра) определяются на основании выражения 


р р 2 
У. -У, 
лови, хр С иг г 
Р;= С — у ‚= К,....0,....К. (22) 
У, —У„, 
ХиРань о.) оф 


Оценка параметров движения космического объекта вычисляется в соответствии с 
формулой 


ОУ ВО (23) 


Оценка элементов вектора интенсивностей маневра осуществляется по формуле 


К 


„=> Ре. (24) 


/— 


Из анализа выражений (15) - (24) следует, что в случае, когда измерения по 
координатам считаются независимыми [3], то рассматриваемый фильтр может быть 
разделен на три независимых фильтра, что приведет к снижению размерности 
обрабатываемых матриц до трех вместо девяти, а это позволит существенно снизить 
вычислительные затраты, необходимые для реализации фильтра. 


Упрощение фильтра 


Рассматриваемый рекуррентный адаптивный фильтр требует для своей 
реализации выполнения достаточно большого количества матричных операций 


вследствие необходимости (+1) раз вычислять о. Оо. иР на каждом шаге 
эп м]? п] п) 


обработки информации о космическом объекте. В работе [3] предложен подход, 
который позволяет упростить адаптивный фильтр. Суть его состоит в том, что 
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находится взвешенная оценка экстраполированных значений параметров, которая 
затем используется в одном обычном фильтре Калмана (12). При этом указанная 
оценка вычисляется по формуле 


а), ВОЕи |, (25) 


Структурная схема такого упрощенного фильтра представлена на рис. 3. 


Ковариаци- 
онная 
матрица 
ошибок 
оценок 
параметров 
=— 


Вычисление текущих коэффициента усиления, 
—> ковариационной матрицы, 
апостериорных вероятностей 


] > 
Оценивание Вектор оценок 
м Блок 
——> координат, скоростей > 
Вектор ИЯ памяти 
измерений и ускорений 
4 
—= 
Экстраполятор (-К) |-> Умножитель -———> 
Нар ЗеВНН 
ЕЙ 
[@ у 
Экстраполятор (0) [= Умножитель -——= чУММатор 
РЕ ЕН 
4 
Гы 
Экстраполятор (К) [^^ Умножитель > 
4 > 
Набор 
фиксиро- 


ванных Набор 


значений апостериорных 
возможных вероятностей 
ускорений 


маневра 


Рисунок 3 — Структурная схема упрощенного адаптивного фильтра 


Результаты моделирования 


Для исследования эффективности применения предлагаемого фильтра осу- 
ществлено статистическое моделирование его работы. Для получения результатов 
произведено 78 реализаций, что обеспечило уровень доверительной вероятности 
97,5 процента. 

Моделирование движения космического аппарата [9] осуществлялось путем 
интегрирования уравнений движения в топоцентрической системе координат методом 
Рунге — Кутта четвертого порядка 

55 28 З . 2 йе 
Х=-ихг -2(20, — то. + ха — то; 


у=-и(у+",)"* —2(х®, —-20,)+(у+",)®: —то,; (26) 


В=- 2? -2(7%, - 0, +20? — то, 
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> а Я 
где геоцентрический гравитационный параметр и = 398603 км / с ; составляющие 
вектора-строки угловой скорости вращения Земли @, = (®.,©,,©,) определяются в 
соответствии с выражениями ©, = ®, с0$ф, зш А,, @, =@,53Шф,, ®, =@, с05ф, с0$ 4, ; 
ф, — широта точки стояния измерительной информационной системы; параметр 


т=х@,+(у+г,)@,+20.; г, — радиус Земли; А, — азимут топоцентрической 


з 
системы координат; угловая скорость вращения Земли @, =0,727211.10* рад/с. 
Для моделирования маневра космического аппарата в уравнения движения 
добавлялись соответствующие компоненты вектора маневра о У.) в необходи- 
мые моменты времени. 
Моделировался маневр космического объекта, осуществляющего движение по 


круговой околоземной орбите радиусом 6771 км. Корректирующий орбитальный 
маневр осуществлялся под действием непрерывной тяги в течение десяти периодов 


обзора радиолокационной станции. Компоненты вектора маневра составили (10 м/ ее. 
10 м/с”, 10 м/с”). 

В адаптивном фильтре произведено разбиение на пять возможных значений 
интенсивностей маневра (-16,—16,—16), (-8,-8,-8)}, (0,0,0), (8,8,8), (16,16,16). 
В качестве характеристики эффективности работы фильтра применялась 
относительная шаровая динамическая ошибка об,,=0,./б„им› При этом 


о, = АХ? +4? +47?. 
Обработка поступающей информации осуществлялась с помощью обычного 
фильтра Калмана, адаптивного фильтра с линейной моделью движения [3], предла- 
гаемого адаптивного фильтра. 
Результаты моделирования показали, что на участке маневра о 


дотн 


не превышала 
значение 0,5 для предлагаемого фильтра, 1,2 — для фильтра [3], 1,9 — для фильтра 
Калмана. При этом фильтр [3] постоянно работал в режиме маневра. На стационарном 
участке обработки измерений величина с 


зн ДЛЯ рассматриваемого фильтра превосхо- 
дила аналогичную величину для фильтра Калмана на 12 процентов, а у фильтра [3] на 


21 процент. 


Выводы 


1. Разработан рекуррентный адаптивный фильтр, который позволяет эффективно 
осуществлять обработку информации о космических объектах на участке орбитального 
маневра. 

2. Особенностью фильтра является применение в составляющих фильтрах или в 
экстраполяторах упрощенного фильтра квадратичной модели движения, которая в 
большей степени соответствует реальному движению космических объектов на 
интервалах их сопровождения измерительными информационными системами. 

3. В основе синтеза лежит принцип разделения. 
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4. Разработанный фильтр показал свое превосходство перед аналогом и обычным 
фильтром Калмана при работе на участке орбитального маневра, обеспечив динамичес- 
кую ошибку на 50 — 70 процентов меньшую, чем фильтры, с которыми производилось 
сравнение. 
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Г.Ю. Грииин 

Метод адаптивно! рекурентно! ф1льтрацй вимрювань параметр!в космчних 06б’ект!в 

У робот! розглянутий синтез адаптивного рекурентного флльтру, який призначений для статистично! 
обробки вимрювань параметрив руху космичних об’екпв на еташ ор@тального маневру у вимрювальних 
1нформашйних системах. В основ! розробки лежить принцип роздлення. У флльтр! врахована особлив1сть 
модел! руху косм!чних об’ектв, що дозволило ефективно здснювати 1дентифкацию параметрив 
траекторй на етати орбтального маневру, а також 1стотно понизити динам1чну помилку супроводження на 
даному еташ. 


1.Уи. Сизт 
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